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Abstract 

In order to extend the recently proposed Monte Carlo Hamiltonian to many-body 
systems, we suggest to concept of a stochastic basis. We apply it to the chain of 
N s = 9 coupled anharmonic oscillators. We compute the spectrum of excited states 
in a finite energy window and thermodynamical observables free energy, average 
energy, entropy and specific heat in a finite temperature window. Comparing the 
results of the Monte Carlo Hamiltonian with standard Lagrangian lattice calcula- 
tions, we find good agreement. However, the Monte Carlo Hamiltonian results show 
less fluctuations under variation of temperature. 



1 Introduction 



Path integral quantization in the Lagrangian formulation and canonical quan- 
tization in the Hamiltonian formulation are two ways to quantize classical 
systems. The Lagrangian formulation is suitable for numerical simulations on 
a computer via Monte Carlo. The enormous success of lattice gauge theory 
over the last two and half decades is due to the fact that the Monte Carlo 
method with importance sampling is an excellent technique to compute high 
dimensional (and even "infinite" dimensional) integrals. 

Unfortunately, using the Lagrangian formulation it is difficult to estimate wave 
functions and the spectrum of excited states. Wave functions in conjunction 
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with the energy spectrum contain more physical information than the en- 
ergy spectrum alone. Although lattice QCD simulations in the Lagrangian 
formulation give good estimates of the hadron masses, one is yet far from a 
comprehensive understanding of hadrons. Let us take as example a new type of 
hardrons made of gluons, the so-called glueballs. Lattice QCD calculations [1] 
predict the mass of the lightest glueball with quantum number J PC = ++ , 
to be 1650 ± lOOMeV. Experimentally, there are at least two candidates: 
/o(1500) and /j(1710). The investigation of the glueball production and de- 
cays can certainly provide additional important information for experimental 
determination of a glueball. Therefore, it is important to be able to compute 
the glueball wave function. 

In the Hamiltonian formulation, one can obtain the ground state energy, but 
also wave functions and the spectrum of excited states. Often, and in partic- 
ular in the case of many-body systems, it is difficult to solve the stationary 
Schrodinger equation. In Ref.[2] we have suggested how to construct an effec- 
tive Hamiltonian via Monte Carlo, which allows to compute the low energy 
spectrum and the corresponding wave functions. The method has been tested 
in quantum mechanics in D=l,2 and 3 dimensions, for the free system, the 
harmonic oscillator and a number of other local potentials. In all cases, the 
exact results were well reproduced [3]. 

If one wants to construct the Monte Carlo Hamiltonian in the case of high 
dimensional systems or for many-body systems using a regular basis like in 
quantum mechanics, one runs into serious troubles. Consider, e.g. a 1-D chain 
of N spin 1/2 particles. For N = 25, the dimension of Hilbert space is d = 
33554432, which is prohibitively large for numerical computations. A physical 
solution of this problem is the construction of a lower dimensional space but 
retaining the important degrees of freedom. This principle is realized, e.g., is 
the folded diagram technique [4]. In recent years, Monte Carlo methods have 
been widely used to solve problems in quantum physics. For example, with 
quantum Monte Carlo there has been improvement in nuclear shell model 
calculations [5]. A proposal to solve the sign problem in Monte Carlo Greens 
function method, useful for spin models has been made by Sorello [6]. Lee et 
al. [7] have suggested a method to diagonalize Hamiltonians, via a random 
search of basis vectors having a large overlap with low-energy eigenstates. In 
contrast to that, in this work we construct a basis, matrix elements of the 
transition amplitude and hence the Hamiltonian via the path integral starting 
from the action. 

We suggest an extension of the idea of the Monte Carlo Hamiltonian with 
the purpose to treat many-body problems: We try to mimic the success of 
the Monte Carlo method with importance sampling when solving Euclidean 
path integrals in the Lagrangian formulation. There one constructs a "small" 
number (say 100 - 1000) of representative (equilibrium) configurations and 
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computes the expectation value of an observable by summing the observable 
over those configurations. In close analogy to those equilibrium configurations 
we suggest here to use stochastically chosen representative basis states (called 
stochastic basis in the following). This will allow us to construct an effective 
Hamiltonian from transition matrix elements between those stochastic basis 
states. The goal is to solve the many-body problem by construction of the 
effective Hamiltonian in such a "model space". One expects the dimension 
of the effective Hamiltonian to be in the same order as the typical number 
of configurations (100-1000) used in computing path integrals in lattice field 
theory. 

We will present results showing that the effective Hamiltonian in conjunction 
with a stochastic basis works in quantum mechanical many-body systems: 
As example we consider a chain of coupled harmonic anharmonic oscillators 
(scalar 4 model), and compute thermodynamical observables. We find that 
the Monte Carlo Hamiltonian gives quite accurate results for the spectrum 
in a finite energy window, and also for thermodynamical observables like free 
energy, average energy, entropy and specific heat in a finite temperature win- 
dow. 



2 Effective Hamiltonian 



Let us briefly recall the meaning of the Monte Carlo Hamiltonian [2]. Using 
Feynman's path integral formulation [8], we consider the transition ampli- 
tude in imaginary time from t — to t — T. Using imaginary time makes 
the path integral mathematically well defined, and renders it amenable to 
numerical simulations by Monte Carlo. Because the effective Hamiltonian is 
time independent, its construction in imaginary time should give the same 
result as in real time. We consider the transition amplitudes for transitions 
between position states. Let {xi, . . . , x N } denote a discrete set of points. Then 
{\xi), . . . , \xn)} forms a basis of position states. We consider the transition 
amplitudes 



where Se denotes the Euclidean action for a given path x(t) going from Xj, 
t = to x h t = T, 



M lj (T) = <x i \e 



HT/h 



Xj > 




(1) 
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S E [x{t)} = S + S V = J dt -mi 2 + V(x) . (2) 

o 

The numerical computation of the matrix elements can be done using 
standard Monte Carlo with importance sampling, by writing it as a ratio of 
two path integrals, thus expressing it as a generalized expectation value of an 
observable O = exp[— S v /f\\. 



M ij {T) = M$>(T) 



f[dx] exp[-S v [x]/h] exp[-S [z]/ft]|^.J 
f[dx] exp[-S [x]/h]\ 3 () 



I &i 



= M§\T) (eMSv[x]/h]) . (3) 

The matrix elements corresponding to the free action So, are known ana- 

lytically. The transition amplitudes M^(T) form a matrix M(T) = [Mij(T)] NxN . 
This matrix M(T) is a positive, Hermitian matrix. It can be factorized into a 
unitary matrix U and a real diagonal matrix D(T), such that 



M(T) = U ] D{T)U. (4) 
Then from Eq.(l), Eq.(4) one can identify 

Ul=<x, t \E^>, D k {T) = e- E l SST '\ (5) 

The k — th eigenvector \E%" > can be identified with the k — th column of 
matrix W. The energy eigenvalues E e ^ are obtained from the logarithm of 
the diagonal matrix elements of D{T). This yields an effective Hamiltonian, 



N 

He f f = Y,\Et ff >Et ff <Eff\. (6) 
fc=i 



3 Regular basis 



Above we have not specified the distribution of the nodes {xi, . . . ,x N }. Sup- 
pose they are equidistantly distributed over the real axis (x i+ i — Xi = Ax = 
const). Also we have been a bit too cavalier in using position states \xi), which 
are not normalizable. To be rigorous, we need to use normalizable states. This 
can be done by constructing "box" -states, which are normalized and localized 
states. Those basis states are denoted by |ej), % = 1, . . . , N. They are defined 
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in position space by (x\ei) = Ax in the interval Li = [xi 
and zero else. Hence we consider the transition amplitudes 



M ij (T) = <e i ,T\e j ,0> 




= Ax / [dx] exp[-S[x\/h\ + 0(Ax 2 ), i,jel,2,...,N. (7) 



Those "localized" and normalized states, corresponding to equidistantly dis- 
tributed nodes, are called a regular basis. 



4 Stochastic basis 

It is evident that the above basis construction becomes prohibitively large if 
we intend to apply this to a system with many degrees of freedom (many-body 
system). For such situations we desire to construct a small basis which gives an 
effective Hamiltonian and reproduces well observables in a low-energy window. 
Why should such a basis exist? The heuristic argument is the Euclidean path 
integral, which, when evaluated via Monte Carlo with importance sampling, 
gives a good answer for the transition amplitude. This is possible by taking 
into account a "small" number of configurations (e.g. in the order of 100 
- 1000). Roughly speaking, the configurations correspond to basis functions. 
Thus we expect that suitably chosen basis functions exist, the number of which 
is in the order of 100 - 1000, which yields a satisfactory effective low energy 
Hamiltonian. Note, however, that this will be the case only when the basis 
functions are chosen in the "right" way. 

How can we construct such a "small" basis? Let us consider first the most 
simple case, i.e., a free particle in D = 1 dimension. Let us take as "large" basis 
the regular basis, described above. Suppose N is large (N » 1). The idea is 
to make a selection guided by the Euclidean quantum mechanical transition 
amplitude. Recall: For the free system it reads 



Note that this function is positive for all x, y, T. It can be used as a probability 
density. We put y = and define a probability density P(x) by 




(8) 
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1 



Geuci{ x , T; 0, 0), 



Z 




dx G E ud(x,T;0,0). 



(9) 



Then we define a selection process as follows: Using a random process with 
probability density P(x) one draws a "small" set of samples {x v \v G 1, . . . , N stoch }, 
where N stoch « N. 

In the case of the free particle, P(x) is a Gaussian, 



In other words, we select {x u } by drawing from a Gaussian distribution. We 
draw N stoc h samples, giving N stoc h states, where N stoc h is considerably smaller 
than N, the original size of the basis. 

Let us give some thought to the question: Is such probability density physically 
reasonable? Firstly, consider the case when T is large. The Boltzmann-Gibbs 
distribution 



projects onto the ground state when T — > oo. For the free system the ground 
state energy is E — 0. I.e., the distribution Pbg(E) has a strong peak at E — 
(when T — > oo). On the other hand, when T — > oo, then a given by Eq.(10) is 
large. Thus the density P(x), from which we draw the broad Gaussian. 

In the limit a — > oo, it becomes a uniform distribution. Now we go over from 
P(x) to P(k), related via Fourier transformation. If P(x) is uniform, then 
P(k) oc 8(k). Thus it gives the energy E k = ^—\k=o = 0, which is the correct 
ground state energy eigenvalue. Thus, in the extreme low-energy regime, the 
distribution P(x) gives a result consistent with the Boltzmann-Gibbs distribu- 
tion. This is a good indication that P(x) will generate an effective Hamiltonian 
useful for the computation of thermodynamical observables. 

Although less relevant for our purpose, it is instructive to look what hap- 
pens in the opposite situation, i.e., when T is small. In the limit T — > 0, the 
Boltzmann-Gibbs distribution is approximately a constant. All energies occur 
with equal probability. Thus a is also small. The distribution P(x) behaves 
like P(x) oc S(x). The Fourier transform yields P (k) = const., i.e. a uniform 
distribution. Then Ek = is distributed like 1/y/E. This is not the same as 
the Boltzmann-Gibbs distribution. But for small T, which means large energy 
E, it is qualitatively comparable to that of Boltzmann-Gibbs. 




(10) 



P BG (E) = -exp[-E T/h] 



(11) 
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Next we ask: What do we do in the case when a local potential is present? 
The definition of the probability density P(x) given by Eq.(9) generalizes 
to include also local potentials. In order to construct a stochastic basis one 
can proceed via the following steps: (i) Compute the Euclidean Green's func- 
tion Ge(x, t; 0, 0), e.g., by solving the diffusion equation and compute P(x). 
(ii) Find an algorithm giving a random variable x distributed according to 
the probability density P(x) and draw samples from this distribution, giving 
nodes, say x v . Finally, one obtains the stochastic basis by constructing the 
corresponding characteristic states from the nodes x v . 

The same goal can be achieved in an elegant and efficient manner via the 
Euclidean path integral. Writing Eq.(9) as path integral yields 



x,T 

J[dy] cxp[-S E [y]/h] 
P(x) = 



0,0 



j^dxj[dy]exp[-S E [y}/h} 



o.o 



Using a Monte Carlo algorithm with importance sampling, like the Metropolis 
algorithm [9], one generates representative paths, which all start at x — 0, 
t — and arrive at some position x at time t = T. Let us denote those paths 
(configurations) by Cj = Xj(t). We denote the endpoint of path Cj at time 
t = T by xf° = Xj(T). Those form the stochastically selected nodes, which 
define the stochastic basis. 

Like we did above for the regular basis, we construct square integrable box 
states, localized around the stochastic nodes. Those states are denoted by 



|ef °), i — 1, . . . , N. In position space they are defined by (x\ef°) = 1/y Axf° 
for x G If° = [xf° — \Axf°,xf + \ Axf°\i and zero else. Hence we consider 
the transition amplitudes 

M lj {T) = <ef°,T\ef°,0> 

x',T 



; / dx' [ dx" [ \dx] exp\—S\x\/h] 

\ r sto\ sto J J 7 L J L L U J 

Y^AXj L-±d,j I?to Isto 



x",0 



3 

to 



[dx] exp[—S[x]/h] 



x sto ,0 



+ 0((Ax sto ) 2 ), i,jel,2,...,N . (13) 



The above expression involves Axf°, i.e. the volume of the interval If°. In 
a 1-dimensional integral the intervals If° have to be chosen such that they 
cover the domain filled by the stochastic nodes. However, in higher dimensions 
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this will be complicated. How should one choose then the volume of such 
"interval" ? This can be answered by recalling how to compute an integral via 
Monte Carlo with importance sampling. 



J = ( dx P{x) g(x) ^Y,Axi P(xi) gfa) . (14) 
a <=i 

Suppose the nodes of integration X; L are drawn from the distribution P(x), 
normalized to unity. Then the Monte Carlo estimator of J is given by 



1 N 



Jest = !TfJ29( X i) ■ ( 15 ) 



This tells us that the volume of the "intervals" is given by 



Ax - = vm ■ (16) 

The result holds in arbitrary dimensions. 



5 The model: chain of coupled anharmonic oscillators 



We consider a one-dimensional chain of N s coupled harmonic oscillators, with 
anharmonic perturbation. Its Euclidean action is given by 



r Ns 1 O 2 O 2 \ 

S = j dt E + \ (0„ +1 - + l ffi + \<& • (17) 
In the continuum formulation it corresponds to the scalar model, 



Introducing a space-time lattice with lattice spacing a s and a t , this action 
becomes 



N s N t -1 
n=l k=0 



1 / $(x n ,t fc+ i) - $(x re ,t fc ) \ 2 

2 V at ~ J 
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+ 



i f $(x n+1 ,t k ) -${x n ,t k y 



m 



+ — $ 2 (x„,t fc ) + j^\x n ,t k ) 



(19) 



The actions given by Eq.(19) and Eq.(17) can be identified by posing <fi = 
^Ja~l<&, = l/a s , Q = m, and A/2 = g/4!. 



5. 1 Estimation of statistical errors 



We have computed the transition matrix elements via Monte Carlo. This yields 
matrix elements with statistical errors SMij. Using stationary perturba- 
tion theory to lowest order, one can compute the propagation of the statistical 
errors into the energy eigenvalues and wave functions of the effective Hamil- 
tonian. Here we have estimated the error propagation numerically. We have 
considered the matrix My ± SM^ and diagonalized it and computed corre- 
spondingly H e ff via Eqs.(4,5). This gives upper bounds on the error of energy 
eigenvalues and on the error in the wave functions. 



5. 2 Spectrum 



In Tab. 1 we present the energy spectrum. 



5.3 Thermodynamical observables 



A solid test of the Monte Carlo Hamiltonian method is a comparison with 
results from standard Lagrangian lattice calculations. However, the strength 
of the latter approach lies not in the computation of excitation spectra. On 
the other hand, it does very well for the computation of thermodynamical 
observables at thermodynamical equilibrium. The information of the energy 
spectrum enters into such thermodynamical functions. Thus we have chosen 
to compute the following thermodynamical observables: the partition function 
Z, free energy F, average energy U, specific heat C, entropy S and pressure 
P. Those are defined by 



Z(P) = Tr[exp(-PH)] , 
F((3) = ~\ogZ , 

U(P) = ±Tr[Hexp(-pH)] = 
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Table 1 

Energy spectrum from MC Hamiltonian. Model parameters: Q = 1, Qo = 2, and 
A = 1 {h = 1, ks = 1). Approximation parameters: N s = 9, a s = 1, (3 = 2. 



n 




statistical error 


n 


pcff 

D Ii 


statistical error 


1 


11.278101 


0.013907 


29 


15.752267 


0.068596 


2 


13.412064 


0.035563 


30 


15.846867 


0.046377 


3 


13.566906 


0.034268 


31 


15.875840 


0.021624 


4 


13.711700 


0.044635 


32 


15.912447 


0.079069 


5 


13.729973 


0.040317 


33 


15.956051 


0.065853 


6 


13.877653 


0.022487 


34 


15.965277 


0.080568 


7 


14.043190 


0.019916 


35 


16.169730 


0.119116 


8 


14.093273 


0.001820 


36 


16.220694 


0.197029 


9 


14.255088 


0.023932 


37 


16.222774 


0.049391 


10 


14.256074 


0.024682 


38 


16.324841 


0.071902 


11 


14.383945 


0.006369 


39 


16.409331 


0.064214 


12 


14.450105 


0.003479 


40 


16.520363 


0.222361 


13 


14.529230 


0.041190 


41 


16.634297 


0.074976 


14 


14.649623 


0.090273 


42 


16.635276 


0.099756 


15 


14.723430 


0.007810 


43 


16.850276 


0.124949 


16 


14.793013 


0.035949 


44 


16.869804 


0.578303 


17 


14.812091 


0.031060 


45 


17.276424 


0.051881 


18 


14.904202 


0.012391 


46 


17.373372 


0.346835 


19 


15.139851 


0.068791 


47 


17.473837 


0.613308 


20 


15.147149 


0.024694 


48 


17.658680 


0.392987 


21 


15.159335 


0.008713 


49 


17.783055 


0.579326 


22 


15.321950 


0.031130 


50 


18.492900 


0.445111 


23 


15.348790 


0.160528 


51 


18.562615 


0.010174 


24 


15.367387 


0.037622 


52 


18.716123 


0.630255 


25 


15.514476 


0.040735 


53 


18.800322 


0.063546 


26 


15.534130 


0.047999 


54 


19.228190 


0.691547 


27 


15.585427 


0.060088 


55 


19.585549 


0.129330 


28 


15.619200 


0.017032 


56 


19.968143 


2.334383 
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S(P) = -(U-F)=k B p(U-F) , 

T 

(20) 

Here k B denotes the Boltzmann constant. The temperature r is related to (3 
via p = T/h= l/(k B r). 

(a) Computation of thermodynamics from the Monte Carlo Hamiltonian: 
When we approximate H by H ef[ , we can express thermodynamical observables 
via the eigenvalues of the effective Hamiltonian 



N 



n=l 



F cS (P) = -^ log Z eff (/5) , 



C eff 03) = £; B /3 2 



1 



S cS ((3) = k B (3 \U CS {(3) - F cii {(3) 



N N 



_n=l n=l 



peS{p) = _ F<S(J3,V + AV)-F<S(J3,V) ^ ^ 

All eigenvalues have been computed from matrix elements, Eq.(7), in which 
enters the transition time T or the corresponding value of (5. When going 
to the continuum limit a s — > 0, a t — > 0, as well as to the thermodynamic 
limit (infinite volume limit), the energy spectrum should become independent 
of the temperature (i.e. of (3 and r) ("perfect scaling"). In practice we have 
worked on a finite lattice volume (N s nodes in spatial direction) and also used 
finite lattice spacing a s and a t . Consequently, the energy spectrum obtained 
will have a finite volume dependence and also some dependence on the lattice 
resolution a s and a t . In practice this manifests itself in the existence of an 
"energy window" and a "temperature window" ( "scaling window" ) where the 
eigenvalues are close to those of the continuum and infinite volume limit, 
and consequently depend very little on (3 entering the transition amplitudes. 
This scaling behavior can be improved and the window can be enlarged by 
approaching the continuum limit, and the infinite volume limit. It also requires 
to increase f3 and to improve the statistics in the Monte Carlo computation of 
matrix elements. The numerical analysis of those scaling properties merits a 
detailed numerical study, which we defer to a future study. In the numerical 



results presented below we have computed the spectrum at a fixed transition 
time T = (3 = 2 (h = k B = 1). 



One should note that once the energy spectrum has been obtained in such a 
window, thermodynamic functions can be computed easily for all values of (5 
in the temperature window, from Eqs.(21). This property is a nice feature of 
the Monte Carlo Hamiltonian approach. It is in contrast to the Lagrangian 
approach, where all thermodynamic functions at a particular value of f3 require 
an independent simulation. 

(b) Computation of thermodynamics in the standard Lagrangian lattice for- 
mulation: 

First, let us consider the average energy U{j3). It is given by the path integral 



,t=o 



bo,t=P 
b ,t=0 



(22) 



Putting f3 = N t a t , differentiation with respect to (3 can be expressed as differ- 
entiation with respect to a t . One obtains 



N t -1 
1 



4> ,t=p 
4> ,t=o 



2a, N, 



IlUto'd^] expi-Stt}] 



<t>o,t=f3 



,t=0 



2a t N t K da t 1 ' 



(23) 



The computation of the free energy F poses some problem in Lagrangian 
lattice formulation. This has to do with the fact that unlike the average energy 
U, F can not easily be written as some expectation value of an observable (like 
O = -^-S[(f>] for average energy). One possibility is to first compute U(/3) and 
then integrate over (3. However, that gives F((3) only up to a constant. Another 
possibility is to write the partition function of the scalar model as partition 
function of the Klein-Gordon model times an expectation value involving the 
c/> 4 interaction in a path integral with the Klein-Gordon action as weight factor, 



Z<t> 4 (P)=Z KG (P) (exph^ 4 ])^ . (24) 
Then one obtains the free energy from 
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F*\(3) 



F KG (P)-^\og(eM-^ 4 ])K G (P). 



(25) 



The entropy, being essentially a difference between average energy and free 
energy is easily obtained from the former two quantities. 

Like the average energy U(/3) also the specific heat C(/3) can be expressed as 
a combination of expectation values involving derivatives of the action with 
respect to a t . 



C{(3) = -k B (3 2 



N. 



2N t a 2 



,d 2 S ( dS_ )2) 



da 2 



(26) 



Let us compare the results from the MC Hamiltonian with those from the 
standard Lagrangian lattice approach. We haven chosen N s = 9 and a s = 1 
(note that we have made no attempt to go to the continuum limit of the 
quantum theory, but our purpose is to compare both methods on given finite 
lattice). As model parameters we have taken Q = 1, Qq = 2, and A = 1 
(h = 1, ks = 1). In the Monte Carlo Hamiltonian simulation we used a 
stochastic basis of N stoch = 100 states, and we used N con f = 300 configurations 
to measure the matrix elements. In the Lagrangian lattice simulation we used 
Nconf = 100 — 200 configurations to measure the free energy F, and N con f = 
100000-200000 configurations to measure U and C. It took about 192 h CPU 
on a Pentium 500 computer to obtain all Monte Carlo Hamiltonian results, 
while the Lagrangian lattice simulations took 2.8 h CPU for each value of (3. 

We found in the Monte Carlo Hamiltonian approach that the matrix elements 
are not very sensitive to the choice of a t , provided that a t is sufficiently small. 
This means the discretization error in the path integral is quite small. The 
results presented in the Figs, below correspond to a t = 1/30. In the Lagrangian 
lattice approach U is stable for this value of a t . However, the results for F and 
C are very sensitive to the choice of a t . Note, in Eq.(25) for the free energy, 
the analytical formula F KG (/3) is has been used, which corresponds to the 
limit a t — > 0. Thus we have chosen a t to be small, a t = 0.01, in this case. On 
the other hand, when measuring the specific heat, a small value of a t leads 
to very large fluctuations. C is stable only for a t > 0.1. Therefore, we have 
chosen a t — 0.1 in that case. 

Fig. 1 shows the free energy F as a function of (3. We make the following 
observations: First, there is good overall agreement in the range 1 < f3 < 10. 
Second, the Lagrangian lattice data fluctuate more than those from the MC 
Hamiltonian. However, the estimated statistical errors are comparable. Fig. 
2 shows the average energy U(f3). The behavior is qualitatively the same as 
for the free energy. One notes that at f3 — 0.5 a marked difference shows up 
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between MC Hamilton and Lagrangian lattice data. Fig. 3 displays the entropy 
S((3). Because the entropy is essentially given by the difference of average and 
free energy, amplified by the factor /3, one observes an amplification in the 
fluctuations of the Lagrangian lattice data when (5 increases (temperature goes 
to zero). For the same reason, also the statistical errors increase. However, 
the MC Hamilton data are stable in that limit. Again a marked difference 
between the two methods become apparent at f3 — 0.5. Finally, we display in 
Fig. 4 the results of the specific heat C(/3). The computation of this function 
in the Lagrangian lattice approach involves second order derivatives and the 
occurrence of cancellations (see Eq.(26)). This requires very high statistics, 
in particular for large j3 {(3 > 3). Hence we have measured C(/3) only up 
to (3 = 3 in the Lagrangian lattice method. From the data of F, U, S and 
C, we estimate the temperature window of the Monte Carlo Hamiltonian to 
range from f3 — 1 to f3 — 10. As discussed above, the size of the window is 
expected to depend on the model parameters as well as on the approximation 
parameters. 



6 Discussion 



We have suggested how to obtain an effective low energy Hamiltonian by 
constructing via Monte Carlo a stochastic basis. We want to stress that the 
advantage of using a stochastic basis shows up in high-dimensional systems 
and in many-body systems. This is based on the experience with integrals, 
where summation over Monte Carlo nodes wins over fixed node rules for di- 
mensions D > 6, as a rule of thumb. We have shown in the scalar model that 
the MC Hamiltonian with stochastic basics works well by computing ther- 
modynamical observables. While thermodynamical functions can be obtained 
also from the standard Lagrangian lattice approach, we found that the results 
from the MC Hamiltonian display less fluctuations. Moreover, we have shown 
that the MC Hamiltonian also provides the spectrum and corresponding wave 
functions in some finite energy window, which is very difficult to obtain in the 
Lagrangian lattice approach. It is here, in our opinion, where the MC Hamil- 
tonian approach has an advantage. Examples of physics where wave functions 
of many-body systems play a role are: Hadron structure functions in particle 
physics, electromagnetic form factors in nuclear physics, Bose-Einstein con- 
densation in atomic physics. We hope that the Monte Carlo Hamiltonian will 
allow to make progress in those areas. 
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Fig. 1. Free energy F{(3). Comparison of results from Monte Carlo Hamiltonian 
(filled circles) with standard Lagrangian lattice calculations (open circles). 
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Fig. 2. Same as Fig. 1, for average energy £/(/?). 
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Fig. 3. Same as Fig. 1, for entropy S(/3). 
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Fig. 4. Same as Fig. 1, for specific heat C{f5). 
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